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Q\ ' Abstract 
O ■ 

Q"^ , A general procedure is presented to construct conditionally exactly solvable 

I (CES) potentials using the techniques of supersymmetric quantum mechanics. The 

• method is illustrated with potentials related to the harmonic oscillator problem. 

Besides recovering known results, new CES potentials are also obtained within the 

framework of this general approach. The conditions under which this method leads 

' to CES potentials are also discussed. 

^3 ■ 



- 

X 



1 Introduction 



^ i Solvable potential problems have have played a dual role since the beginnings of quantum 

mechanics. First, they represented useful aids in modelling realistic physical problems, 
and second, they offered an interesting field of investigation in their own right. Related 
to this latter area, the concept of solvability has changed to some extent in recent years. 
Besides exactly solvable problems, for which the bound-state energy spectrum and solu- 
tions could be given in general analytical form, quasi-exactly solvable (QES) (see e.g. 
and conditionally exactly solvable (CES) potential classes have also been identified 
recently. 

In the first case only a finite number of eigenstates can be obtained exactly, while 
in the latter one analytical solutions are available only if some (or all) of the potential 
parameters are fine tuned to specific numerical values. 
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There are different types of CES potentials, which is also reflected in the way they 
can be most naturally constructed. Some of them, including the first CES potentials 
P, H (see also ref |Q for some interesting comments on the construction of CES potentials) 
have their bound-state solutions in terms of a single special function (polynomial); a 
structure characteristic of Natanzon class potentials For another class, these solutions 
have more complex structure, but generally they can be expressed in terms of the linear 
combination of two special functions. Typically these CES potentials are constructed as 
supersymmetric partners of some simple potentials 0]. Their CES nature hinges on 
the fact whether the parameters of their partners can be chosen such that they can be 
reduced to some simple potential with known solutions and energy eigenvalues. According 
to the techniques of supersymmetric quantum mechanics, or SUSYQM (for reviews see e.g. 
1^), the CES potentials constructed in this way are then essentially isospectral with their 
partners, i.e. the two spectra are identical or differ only in their ground state. The bound- 
state solutions of CES potentials are obtained from those of their simple (Natanzon-type) 
partner potentials by acting on these latter ones with linear differential operators. This 
explains why their bound-state solutions possess the specific structure described above. 

In Refs. 1^, 1^ some CES potentials have been constructed by SUSYQM. The aim of this 
publication is to show that this procedure can be made more systematic by making use of 
various types of SUSYQM transformations. The rather general nature of this treatment 
allows the recovery of known results and also the derivation of new CES potentials in 
the same framework. Our examples concern CES potentials related to the harmonic 
oscillator potential in three or one dimension (the standard examples of Refs. 0]), but 
the formalism is equally applicable to other types of potentials as well. 



2 The conventional SUSYQM approach to CES po- 
tentials 



Let us assume that there is a pair of SUSYQM partner potentials V± (r), which can be 
constructed from a superpotential Wo{r) in the usual way: 



Vi'\r) = l[W^{r)±W;,{r)]. 



(1) 



Consider now a superpotential of the form 



W{r) = Wo{r) + w{r) . 



(2) 



The partner potentials generated from W{r) are then 




(3) 
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V_(r) = Vi^\r) + Wo{r)w{r) + ^[w\r) - w'{r)] . (4) 

Let us now insist on that one of these potentials, say V-^.{r) is related to some known 
potential up to an energy shift. In the simplest case this could be vf'\r) in Eq. (P: 

V+{r) = vi''\r)+A . (5) 

Combined with this requirement immediately introduces a Riccati-type differential 
equation for w{r): 

^[w^{r) +w'{r)] + Woir)w{r) = A . (6) 

If this equation is solved, then a pair of SUSYQM potentials is obtained, from which one 
of the partner potentials, V+{r), corresponds to a known potential (up to an energy shift). 
Therefore, both the spectrum and the wavefunctions of the partner potential V^{r) can 
be obtained in the usual way. 

In the examples in ^ vf^\r) was the harmonic oscillator potential in 1 and 3 dimen- 
sions, with Wo{r) being the corresponding superpotential. In both cases the structure of 
w{r) was of the type 

In the practical examples N=l was used. The difference was that in the one-dimensional 
case the authors of [^] considered unbroken supersymmetry (and therefore V-{x) had one 
more state than V+(x)), while in the three-dimensional case they chose to discuss broken 
supersymmetry, so the spectra of the partner potentials were identical. 



3 An alternative SUSYQM construction of CES po- 
tentials 

Here we use various SUSYQM transformations systematically to recapitulate the formal- 
ism of Section ^ and to put it into a more general context. 

It is known from the theory of isospectral potentials that a potential V2(r) isospectral 
with a known potential Vi(r) can be constructed by 

V^2(r) = Vl(r)-^ln0(r) , (8) 
where 0(r) is a solution of the Schrodinger equation 

-l^ + V^ir^ = ^^- (9) 
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Here e is usually called factorization energy. Depending on the value of e and the boundary 
conditions of the solution 0(r), V2(r) in (^) will have various properties. In the case of a 
radial problem (in three-dimensions) four types of transformations are possible. These are 
related four different types of nodeless solutions 0(r) of Eq. and have been described 
in terms of SUSYQM |10|, |12|. The nodelessness of 0(r) guarantees that the result- 



ing potential V2(r) does not have singularities for finite values of r (besides the origin), 
and this can be achieved whenever the factorization energy e is below the ground-state 
energy of Vi{r) [^j. We briefly summarize the basic characteristics of the four SUSYQM 
transformation types in Table |l[ 

Let us consider the radial harmonic oscillator as an example and solve (^ for 0(r) with 

Vi(r) = l/f(r) = l{WS(r) + Wi(r)] 



1,^ + 2(1+1) +, + 5. (10) 



Here the superpotential is Wo{r) = r + (7 + l)r ^, and the bound states of Vi(r) are found 
at En = 2n + 2'y + 3. The solution 0(r) can be searched for in the form 

(f){r) ~ exp F{a, b; Cr'^) , (11) 



where F(a, b; z) is the confluent hypergeometric function ||T3|. Straightforward calculation 
shows that (^ transforms into the confluent hypergeometric equation if the following 
conditions hold: 

A(A-l) = 7(7 + 1) , 52 = 1, B = -C- (12) 
b = A+^- (13) 
e 7 3 A 1 

Recalling that besides F{a, b; z), z^~^F{a — b + l, 2 — 6; z) is a linearly independent solution 



of the same confluent hypergeometric function |]T3[, the general solution (j)(r) has the form 



Hr) . exp (f ■^) I a,r--^FiJ^ + ^ + 1 + 1,^1 Cr^) 

Note that the two terms in Eq. (|1^) are connected by the 7^1—7 transformation, 
therefore it is enough to consider one of the solutions (A = 7 + I or A = —7) of A{A — 1) = 
7(7 + 1) in Eq. (0). The solutions corresponding to the transformations Ti, T2, T3 and 
T4 in Table |l] can then be identifled by imposing the appropriate boundary conditions on 
0(r). 
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Substituting the 0(r) function in Eq. (|^) one obtains an expression for V2(r) in terms of 
0'(r) and 0"(r). With the use of V2(r) can be expressed in terms of Vi(r), e and 
In this last expression the first-order derivatives of two confluent hypergeometric func- 
tions occur, each of which can be expressed in terms of another confluent hypergeometric 
function |]13[. This means that V2(r) can be expressed in a somewhat complicated, but 
closed analytic form. A special situation occurs when a = —N or a — 6 + 1 = —M holds. 
In this case one of the confluent hypergeometric functions occuring in ([T5|) reduces to an 
iV-th or M-th order (generalized Laguerre |jl3|) polynomial of the argument. According 
to (p!4|), this case corresponds to specific choices of the factorization energy e. We note 
that in principle both confluent hypergeometric functions can reduce to a polynomial if 
both e and 7 has a specific values: e = — C(M + A^ + 1)— 7 — |, 7 = M — iV-|-|. Let us 
now consider the four transformations Ti, T2, T3 and T4 one by one. 



3.1 The Ti case 

The boundary condition at r = allows solutions only with q;2 = in ([T5|) (if 7 > 
holds), furthermore i? = — 1 is also required to fulfil the asymptotic boundary condition. 
The a = —N condition leads to e = 2N + 27 + 3. This factorization energy corresponds 
to the bound-states energies of Vi(r) = V^\r) and (j){r) simply reproduces the physical 
wavef unctions. It is known that for N = n ^ the transformed potential V2(r) has 
singularities at those locations, where the wavefunctions have nodes. The ground-state 
wavefunction with n = 0, however, is nodeless, and the Ti transformation then simply 
retrieves the classic SUSYQM transformation which eliminates the ground state of Vi(r) 
and increases the value of 7 with one unit. 



3.2 The Tg case 

Similarly to the Ti case only the term regular at the origin is allowed by the boundary 
condition at r = (i.e. 02 = 0), however, the asymptotic boundary condition requires B = 
1 in this case. The a = —N polynomial condition then leads to the specific factorization 
energies e = —2N, which are always below the ground-state energy of Vi(r), so the 
nodelessness of (f){r) is always secured. The = choice recovers V2(r) as another 
oscillator with the same spectrum as Vi{r): only the value of 7 is increased with one unit 
and the energy is shifted downwards with one unit. The A^ = 1 case results in the CES 
potential described in (denoted by V_{r) there) up to an energy shift: 



1 2 , (7 + l)(7 + 2) , , 1 , 2g, . 
V2{r) = -r + — +7 + 7^ + 77— ^-TT 2' ^1 = I Q • U6 

2 2r^ 2 (1 + Qir-^Y 1 + gir-^ 27 + 3 



The energy shift is two units here, and it appears both in the numerical constant in ( [I6| ) 
(it is 1/2 here and 7/2 in 0) and the factorization energy (—2 here and in P]). Similar, 
but more complicated isospectral potentials would arise from choosing A^ > 1. 
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3.3 The case 



In contrast with the previous two cases, the boundary condition at the origin now allows 
both the regular and the singular solution in (|T5|). The ratio of the two coefficients, ai 
and a2 appears as a new parameter in V2{r): the resulting potential family will have the 
same spectrum, but different shape. Similarly to the T3 case, this one is usually also 
interpreted as a situation with broken supersymmetry, because the spectra of the partner 
potentials (and, of course, of the whole family) is identical. 

In order to get a situation similar to the T3 and Ti cases, we can restrict this potential 



family to a single potential, i.e. to that with ai = in (]T5|). The structure of 0(r) 
then becomes the same as before: it will contain only a single confluent hypergeometric 
function. With the loss of the generality, however, the nodelessness of 0(r) cannot be 
guaranteed in general, rather it has to be checked in each case separately. 

Taking also into account the asymptotic boundary condition which now requires B = 
— 1, we find that the a — b + 1 = —N polynomial condition now leads to factorization 
energies e = 2N + 2. The = choice again results in another harmonic oscillator 
potential, with 7 decreased with one unit and with an energy shift of one unit upwards. 
For N = 1 a. potential similar to that in P] arises, whenever 7 > 1/2 holds. (As we have 
mentioned already, this latter condition secures that the polynomial F{—1, —7 + |; r^) = 
1 + 2r^/ (27 — 1) remains nodeless, and there will be no singularities in the V2(r). In fact, 
the 2'y + 3 = Eq > e = 2N + 2 condition also leads to 7 > 1/2 for = 1.) The potential 
is then 



T./ N 2 7(7-1) 5 45(?r2 2gi 



„2 

2 ' 2r2 ' ' ' 2' (l+^ir2)2 I + g^r^ ' 27 - 1 



The functional form of V2(r) is essentially the same as that of (p!6|), only the value of 7, 
the numerical constant and gi is different. Similarly to the T3 case, further potentials 
isospectral with a harmonic oscillator can be constructed by choosing A^ > 1, but the 
nodelessness of (hir) has to be checked in each case. 



3.4 The T2 case 

The situation here is the same as in the T4 case: both the regular and the singular solutions 
are allowed by the boundary condition at the origin. This means, that we again have a 
whole family of potentials V2(r), which have the same spectrum and differ only in their 
shape. As before, we again restrict our attention to a particular member of this family, 
i.e. to the potential obtained with ai = 0. Furthermore, we consider the polynomial 
condition a — b + 1 = —N, which leads to e = — 2A^ + 27 + 1. 

For A^ = 0, V2{r) = ^r^ + "'^Jr^^^ ~^ I' which corresponds to another harmonic oscillator 
potential with the 7 value decreased by one unit and also shifted lower by one energy 
unit. Clearly, this corresponds to the usual SUSYQM transformation which inserts a 
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new state (at E = 2'y + 1) below the ground state of Vi(r). For = 1 we find that 
F{ — 1, —7 + — = 1 — 2r^/ (27 — 1), which has a node at a positive value of r, unless 
7 < 1/2 holds. As in the T4 case, here we have to check the nodelessness of (f){r) in 
each case, because it cannot be automatically guaranteed after we restricted the general 
solution by selecting ai = in (pTSf ). The functional form of V2(r) is the same as (p!?!) , but 
with gi = 2/(1 - 27). 



3.5 The one-dimensional case 

We note that similarly to the radial equation, the one-dimesional case can also be handled 
in the present framework. The difference arises from the different boundary conditions. 
The various transformation types corresponding to solutions of a Schrodinger equation 



with different boundary conditions has been described in [|rT[|. Here we only mention the 
example discussed in for one dimension. In order to construct V2(r) with one more 
bound state than Vi{r), the solution diverging in both directions has to be considered 



11]. In general, such a solution should be constructed from the two linearly independent 



solutions of the Schrodinger equation, similarly to the situation seen in Subsections and 



3^ . However, taking only one of these, 0(a;) = exp(^)F(— 1, |; —x"^) = exp(^)(l + 2x^ 

we obtain the new state introduced for V2{x) as ^^o\x) ^ l/</)(x). Note that V2 (r) is 
symmetric with respect to the x ^ —x transformation. In the general case V2(r) would be 
asymmetric, similarly to the SUSYQM partner potentials of the one-dimensional harmonic 



oscillator constructed in Ref. M . 



4 Discussion 

The relation of the two procedures outlined in Sections ^ and |^ can be interpreted in a 
simple way by noting that the partner potentials are linked by V^{r) — V^{r) = W'{r) 
and Vi{r) — V2{r) = {ln(f){r))". From this 

W{r) = {\n(f){r)y + c (18) 

follows. Direct integration of @ and (0) with Wo{r) = r + (7 + l)r~^, as in and c = 0, 
indeed, recovers the general solution 0(r) specific to the T3 case: 

0(r) ^ r^+i exp (^^^ Ul,{l + g,r') . (19) 

In addition to the notation of , (70 = was also introduced for convenience. This function 
is also an A^'th order polynomial, as expected from ([TT| ) for a = —N. In addition to the T3 
case, the situation should be the same for the other cases mentioned here, including also 
the one-dimensional case. We note that in the T2, T4 and in the one-dimensional cases 
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discussed here, (|19D is not the most general form of the solution, rather it is a specific 
member of a family of solutions obtained as the linear combination of two independent 
solutions. 

In fact, all the V2(r) potentials derived from polynomial-type 0(r) solutions can be 
expressed by a common formula. Substituting a = —N in (1TT|) and combining it with (|^) 
and (0) one gets 

The solutions relevant to the Ti, T2, T3 and T4 cases can then be obtained by substituting 
[A, B, C] = [7 + 1, -1, 1], [-7, 1, -1], [7 + 1, 1, -1] and [-7, -1, 1], respectively. In the 

= case the last term in ( ^0|) cancels and V2(r) contains only terms characteristic of the 
three-dimensional harmonic oscillator potential. For = 1, 1, ^ + |; Cr^) = l + gir"^, 
with Qi = —2C/{2A + 1), which gives rise to two new terms, formally identical ith the 
last two terms of (p^G]). 

Another interesting formula can also be derived if recalling (|^) and (p!0|), i.e. Vi(r) = 
K|_(r) — A, which also indicates V2(r) = V_{r) — A. From this Vi{r) + V2{r) = W'^{r) — 2A 
directly follows. On the other hand, (H) and (|) lead to V2{r) + Vi{r) = 2e + (070)^. 
Combinig these two formulas, we get 

[W{r) - (In 0(r))'] [W{r) + (In 0(r))'] = 2(e + A) . (21) 

This means, that W{r) = ±(ln(/)(r))' always requires e = —A to hold. In this way A, the 
constant appearing in the Riccati equation (^ in Section |^ is related to the factorization 
energy used in the formulas in Section ^. Note that the W = — (ln0)' choice is also 
acceptable here, and it simply corresponds to the V+ ^ V+ and Vi ^ V2 replacements. 

We stress that although the examples we presented here are related to the harmonic 
oscillator in three and one dimensions, our treatment can be applied to other types of 
potentials as well. Besides the relatively simple shape-invariant potential [|15[ (see also 



iTq , for a list of them), any member of the more general Natanzon potential class [g can 
also serve as the V^\r) reference potential in (^). However, for the first such applications 
some more thoroughly studied special Natanzon potentials could be the best candidates 
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limr_>oo0 convergent 
Spectrum deletes ground s 

modification 

Singularity (7 + l)r~^ 

modification 



T2 


T3 


T4 


e < Eq 


e < Eq 


e < Eq 




^7+1 


r-T 


divergent 


divergent 


convergent 


new ground state 


none 


none 


(0 < 7 only) 

_2 

— 7r 


(7 + l)r-2 


(0 < 7 only) 
— 7r ^ 
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